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Abstract. We extend the method of Ghasemi and Marshall [SIAM. J. Opt. 
22(2) (2012), pp 460-473], to obtain a lower bound / gPj M for a multivariate 
polynomial /(x) g R[x] of degree < 2d in n variables x = (xi, . . . , x n ) on the 
closed ball {x £ M." : ^Zxf d < M}, computable by geometric programming, 
for any real M . We compare this bound with the (global) lower bound / gp ob- 
tained by Ghasemi and Marshall, and also with the hierarchy of lower bounds, 
computable by scmidcfinitc programming, obtained by Lasscrrc [SIAM J. Opt. 
11(3) (2001) pp 796-816]. Our computations show that the bound / gp ,M im- 
proves on the bound / gp and that the computation of f gP} M, like that of / gp , 
can be carried out quickly and easily for polynomials having of large number 
of variables and/or large degree, assuming a reasonable sparsity of coefficients, 
cases where the corresponding computation using semidefinite programming 
breaks down. 



1. Introduction 

Computing a lower bound on the global minimum on W 1 of a multivariate poly- 
nomial is a standard problem of optimization with many potential applications. 
In the last decade, results in polynomial optimization combined with scmidchnite 
programming (for sums of squares representation), have permitted to make some 
progress. For instance, one may compute a lower bound of polynomial / e R[x] on 
E": 

• by solving the problem / sos := sup{A : / — A is sos}, which is a single 
semidefinite program 

• by applying the hierarchy of semidefinite relaxations to the polynomial 
optimization problem inf{/(x) : V/(x) = 0} (assuming that the infimum 
is attained) 

• by applying the hierarchy of semidefinite relaxations to the polynomial opti- 
mization problem inf{/(x) : ||x|| 2 < M}, for sufficiently large M (assuming 
that a global minimum satisfies that bound constraint). 

All those approaches are very powerful and provide good bounds and sometimes 
the exact value. However, so far, and in view of the present status of semidefinite 
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programming, those methods are limited to small to medium size problems, ex- 
cept if some structured sparsity is present (in which case specialized semidefinite 
relaxations can be implemented; see e.g. [7]). 

This limitation of semidefinite programming to implement sums of squares (SOS) 
representations, was the motivation for providing other SOS certificates and yielded 
the sufficient conditions of [4] and subsequently of [2, 5]. And in a recent work 
Ghasemi and Marshall [G] have shown how to compute a lower bound on the global 
optimum of a multivariate polynomial on R™, by solving a certain geometric pro- 
gram. This formulation as a geometric program is based on the sufficient condition 
for a polynomial to be a sum of squares given in [5] , which generalizes the sufficient 
conditions of [2, 4]. Geometric programming (GP) is a convex optimization prob- 
lem that can be solved efficiently for relatively large scale problems. In Boyd et 
al. [ ] it is claimed that GP problems with up to 10 3 variables and 10 4 constraints 
can be solved via standard interior point methods. For sparse GP problems, i.e., 
where each constraint depends only on a small number of variables, the size limit 
can grow up to 10 4 variables and 10 6 constraints! So the interest of the geometric 
programming formulation is that one may now handle polynomials with a large 
number of variables and high degree, especially when the support of / (i.e., the set 
of non zero coefficients) is small. 

Contribution. Our contribution is to extend the geometric programming for- 
mulation of Ghasemi and Marshall [ ] to provide a lower bound on /*,jvf := min{/(x) : 
^2i%i d < M}. The latter problem has its own interest and also serves as an auxil- 
iary problem to provide a lower bound on /* := min{/(x) : x £ R™} when a global 
minimizer is "guessed" to belong to the ball {x : J2i x 1 < M}. Again, and as 
for [6], the main interest of this approach is to be able to handle polynomials with 
large number of variables and/or large degree for which so far, there is no such 
algorithm. Notice that even for a small number of variables, the SOS approaches 
cannot handle polynomials with large degree. 

2. Main result 

Notation and definitions. Let K[x] be the ring of polynomials in the variables 
x = (sci, . . . , x n ), and for d £ N, let M[x]d C M[x] be the vector space of polynomials 
of degree at most d. Let NJJ := {a £ N™ : |a| < d} where \a\ = J2i a i f° r every 
a £ N n . 

Assume now that d > 1. Let e, := (6a, • • • , 6i n ) £ N n , with Sij = 1 if i = j and 
otherwise, and given f = Y1 /a x " € R[ x ]2ti, let: 

Sl(f) := {a £ N? d : f a ? 0} \ {0, 2de x , ■■■ , 2de n } 
A(f) := {a £ Q(f) : f a x a is not a square in M[x]} 
A(/)< 2d := {a £ A(f) : |a| < 2d}. 

Denote the coefficient fid^ by fis.,i for i = 1, . . . , n. 

We first recall the following result of Ghasemi and Marshall [6]. 
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Proposition 2.1. ([(>, Corollary 3.6]) Let f G M[x] 2( 2 and let p be the optimal value 
of the program: 



(1) 



p = min (2d — \a\) 

q£A(/)< m 



2d 



l/(2d-M) 



s.t. y ' z Ui i < / 2di , i — 1, 

«6A(/) 
2d 



(If) (if) =1; aeA(f),\a\=2d. 



where for every a G A(/), the unknowns z a — (z a .i) G R" satisfy z a ^ = if and 
only ifcti = 0. ffere, := n?=i and (^) Q := EKLi toftft ffte 

convention 0° = 1. T/ien /(x) > /(0) - p /or all x € R™. 

The most interesting case is when > 0, z = in which case the 

program (1) is a geometric program. Somewhat more generally, if V i = 1, ...,n 
either (fids > 0) or (fid,i = and on = V a € A(/)), then the program (1) is 
a geometric program. In the remaining cases the program (1) is not a geometric 
program, the feasibility set of (1) is empty, and the output p is oo. 

Problem statement. Let / G R[x] and, for M > 0, consider the problem: 

n 

(2) P M : f*,M ■= min {/(x) : £>f < M}. 

i=l 

Problem Pm has its own interest but is also an auxiliary problem for the uncon- 
strained problem P^ : /* = min{/(x) : x G R n }, when a global minimizer is 
guessed to belong to the ball Bm '■= {x : Y^i x i — M}. Also, notice that the 
sequence (/*,m), M g N, provides a monotone nonincreasing sequence of upper 
bounds on /„ that converges to /* in finitely many steps whenever P^ has an 
optimal solution x* G R n . 

Main result. With M > fixed, to compute a lower bound for /*,m, let A > 
and consider the polynomial fx G R[x] 



(3) 



x i ^ / A (x) = /(x) - A (A/ - £ A > 0. 



i=i 



Lemma 2.2. Lef / G R[x], deg / < 2d and let f x G R[x] be as in (3). Then: 



(4) 



f* jvj > max min f\fx) = max G(A). 

" A>0 x€l" A>0 



G(A) 



Moreover, if either /„ = /*.m or / ^ convex then equality holds. 
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The proof is standard and will be omitted. Actually, one can show that 



max min / A (x) = min/ Al (x) — {fxi )*,M 

\>0 x£l" x€l" 

where Ai is the least A > such that f\ achieves its global minimum on the ball 
B M . Obviously, /*, M > (fxJ^M- If /* = f*.M then 71 = and /* jM = (/aJ*,m- 
If / is convex then / 7 is convex for each 7 > 0. If / is convex and 71 > then 
the minimum of / 7 on Bm is achieved on the boundary of Bm for < 7 < 71, so 
f*,M = (/ai)*,m holds in this case too. 
Note that equality in (4) fails in general. 

Example 2.3. Let n = 1, 2d = 4, f(x) = 2x 2 (x-2) 2 + (l-x 4 ) = x 4 -8x 3 + 8x 2 + l, 
M = 1. Then f„ m = 1 and Ai = 1 so max min f\(x) = min fx, (x) = 0. 

A>0 x£R x£R' 

Observe that for every A > 0, 
(5) G(A) = min/ A (x), 

and so if for every A > 0, G(A) is a lower bound on G(A), then 



(0) 



f* ai > max G(X) > max G(X). 
\>o \>o 



After relabeling if necessary, we may and will assume that 

J2d,l > J2d,2 > ' ' ' > J2d,n- 

The main result of our paper is as follows: 

Theorem 2.4. Let f G E[x], deg/ < 2d. Then 

f*,M > /(0) + M/ad.1 - pjif, 
twi/j feeing f/ie optimal value of the geometric program: 

2d 



(J) 



Pm = min Mtxi + > (2d — 

Zcv,U ^ ' 

a£A(/)< M 



l/(2d-|a|) 



S.t. 



a S A(/) 
/ \ 2d 



< 1, 



1, 



1; a e A(/), |a| = 2d. 



(*) 



(**) 



f2d, 



i < 1 



f2d,i~l—f2d,i 



< 1, i = 2,...,ra, 



and where for every a € A(/), </ie unknowns z a = (z a> i) £ R" satisfy z a ,i = if 
and only if ai =0. 
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A detailed proof can be found in §5. Observe that the difference between the 
programs (1) and (7) is the presence of the constraints (*) — (**) in the latter, which 
reflects the new contribution of the monomial terms Xxj d in the polynomial fx. 

The geometric program (7) is not a direct application of Proposition 2.1 to the 
polynomial f\ to obtain a lower bound G(X) on G(A), followed by a maximization 
with respect to A. Indeed, this leads to the constraint (**) in equality (instead 
of inequality) form, and so (7) would not be a geometric program; however, in 
the proof we show that this equality constraint can be relaxed to an inequality 
constraint as in (7). 

3. Comparison with other bounds 

Comparison with bound of Ghasemi and Marshall. Assume that / € K[x]2d, 
d > 1. As in [(>] we define f gp to be / gp := f(0)—p, the lower bound for /* obtained 
in Proposition 2.1. We also define f gp .M to be / gp ,M := /(0) + Mf 2d ,i — Pm, the 
lower bound for /*,m obtained in Theorem 2.4. Note that the feasible set of (7) 
is nonempty (i.e., / gp ,M is a rea l number), whereas the feasible set of (1) may be 
empty (i.e., / gp = — oo), even in the case where each f 2 d,i is strictly positive. 

Proposition 3.1. 

(1) /gp,M > /gp- 

(2) If M' < M then f gp . M , > f gp>M . 

( 3 ) /gp = J™ /g P ,M- 

Proof. (1) If the program (1) in Proposition 2.1 has no feasible solutions then 
/ gp = — oo so /gp,M > / gP - Suppose now that (1) has a feasible solution z. In 
particular, f 2d ^ > for i — 1, . . . , n. Fix 6 > 0. Then (z, u) with ut = f 2 d.i + $ 
for all i = 1, . . . ,n, is feasible for the program (7) in Theorem 2.4. This implies 
Pm < M(f 2 d, 1 +S) + p for all 6 > so p M < Mf 2dA + p and 

/gp,M = /(0) + M/ 2 d,i - PAf 

> /(0) + M/ M ,i - Mf 2dA -p = /(0) - p = / gp . 

(2) Suppose M' < M. Observe that the set of feasible solutions for (7) does 
not depend on M. Let (z,u) be a feasible solution of (7). Since M' < M and 
u i > /2d,i it follows that M'(ui — f 2 di) < M(ui — f 2d .i)- This implies that 
p M ' - M'f 2 d,i < Pm - Mf 2 d,i, so 

fgptf' = /(0) + M'f 2dtl - p AP > /(0) + Mf 2d ,x - PM - /gp,M. 

(3) It remains to show that if there exists a real number N such that pm — 
Mf 2d .i < N for each real M > then p < N. Suppose p M - Mf 2d ,i < N for all 
M > 0. Then for each real e > there exists a feasible solution (z, u) = (zm, um) 
of (7) such that 

(8) M(«i-/ M|1 )+ (2^ - |«|) 

aeA(f)< 2d 



2d 



<N + e, 
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for M = 1, 2, • • • . As explained in the proof of Theorem 2.4, we may assume 
Ui-f2ds = Uj-f2d,j for all i, j = 1, . . . ,n. Let A = \ M = u 1 -f 2 d,i, so m = / 2 d,i + A 
for i = 1, . . . , n. From inequality (8) we see that MX < N + e, so A — > as M — » oo. 
Since < z ai j < = / 2t 2,i + A, the sequence (z, u) = (z M ,u M ) is bounded so it 
has some convergent subsequence converging to some (z*,u*). If z* is a feasible 
point of the program (1) then we see by continuity that 

l/(2d-H) 

< N + e 



E 

o6A(/)< M 



^ (2d-H) 



j. \ 2d / 



2d J V Z 



so p < N + e and we are done. The fact that z* is a feasible point for (1) is more 
or less clear, by continuity, except possibly for the fact that Ui > =>■ z* Q > 0. 

If \a\ = 2d this follows from the equation f |=J (^) = 1 which, since the z a) , 
are bounded, implies that the z a .i such that a; > are bounded away from zero. 
Similarly for \a\ < 2d the inequality (8) implies that the z a ^ such that cii > are 
bounded away from zero. □ 

Comparison with bounds of Lasserre. Recall that 



/ sos :=sup{A:/-Aey>[x] 2 } 



The inequality /* > / sos is trivial. The inequality / sos > / gp is established in [6, 

Corollary 3.6]. As explained in [3], / sos is computable by semidefinite programming. 

(k) 

Similarly, for each real M > and each integer k > define / sos M to be the 
supremum of all real numbers A such that 

f-\ = a + r(M-Y / ^ d ) 

for some a, r g ^]R[x] 2 , deg(cr) < 2k + 2d, deg(r) < 2k. As explained in [3], the 
sequence ff££ M' k — 0,1, ■ ■ ■ is nondecreasing and converges to /*.m as k — >• oo 
and each f^} M is computable by semidefinite programming. 

Proposition 3.2. f^ } s M > f gp , M - 

Proof. By the proof of Theorem 2.4, / gP) M = maxG(A) where G(A) := {f\) gp . By 
[ , Corollary 3.6], (f\) B os > (A)gp- Thus for any real e > there exists A > 
such that (/a)sos > (fx)sp — fgp,M ~ e ; so t nere exists a € J^IR[x] 2 such that 
/a - (/ g p,M - 2e) = a, i.e., / - (/ gp , M - 2e) = <r + A(M - E A d )- It follows that 
/soi.M - /gp,Af - 2e - Since e > is arbitrary it follows that f^ s M > f sp<M - □ 

Remark 3.3. 

(1) According to [ , Cor. 3.4], |0(/)| = 1 =>■ / gp = / sos = The same is true 
(trivially) if |f2(/)| = 0. Thus if |f2(/)| < 1 and / achieves its global minimum in 
the ball P>m then 

/* = f*,M ^ /sos.M — fgP,M ^ /gp = /sos = /*) 
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SO 



f*,M — fsos.M ~ fgP,M ~~ fgP — /sos — /*• 

(2) There are explicit formulas for / gp and / gp ,M if |A(/)| = 0. Suppose |A(/)| 
0. As usual, we suppose that /2d,i > • • • > /2<2,n- Then 



/gp — 



/(0) if/2d,n>0 
-00 if / 2( j,n < 



and 

7(0) if/ 2d ,„>0 



/gp,Af 



f(0)+M/ 2d , n if/ad,„<Q 



(3) There are also explicit formulas for / gp and f gp .M if |A(/)| = 1 and fid.i = 1, 
i = 1, ... ,7i. Suppose that |A(/)| = 1 and ji6,,i = 1, i = 1, ■ ■ ■ ,n. Let A(/) = {a}. 
There are two cases to consider: 

Case (i). Suppose \a\ = 2d. In this case 

f/(0) if(!t) 2 V»<i 

SP \-oo if (^) 2d a a > 1 ' 

and 

7(0) if(|-) 2 V><l 
J(0) - M ■ ([(^) M a«] 1 / 3U - 1) if (||) 2d a« > 1 ' 
Case (ii). Suppose |a| < 2d. In this case 

/gp = /(0) - [2d - H][(^) 2 V*] 1/(2cHq|) , 

/(0) - [2d- |a|][(^) M a a ] 1 /(M-H) if M > |a| • [(||)2<i Q ,a]i/(2d-|a|) 
/(0) + M- |/a|[(jf|) |Q| a Q ] 1/2d if M < |a| • [(^) M a a ] 1 /(M-|a|) ' 



/gP,Af 



and 

/gp,M 



Example 3.4. Suppose n = 1, 2d = 6, / = a; 6 + 3x 4 — 9x 2 . Applying Remark 
3.3(3), Case (ii), we see that / gp = -2 • 3 3 / 2 « -10.3923 and 



/gP,Af 

In this example one checks that /* = —5, and 

f*,M = 



-2 • 3 3/2 if M > 3 3/2 
M - 9M 1 / 3 if M < 3 3 / 2 



-5 if M > 1 

M + 3Af 2 / 3 - 9M 1 / 3 if M < 1 
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4. Numerical computations 

To compare the running time efficiency of computation of / gP ,M using geometric 
programming with computation of fj® M using semidefinite programming, we set 
up a test over 10 polynomials for each case to keep track of the running times. 
The polynomials considered had highest degree part ^ x\ d with the lower degree 
coefficients randomly chosen integers between —10 and 10, and M was taken to 
be a random integer between 1 and 10 5 (Table 1) . The source code of the Sage 
program to compute / gp ,M an d f^l M , developed by the first author, is available at 
http://goo.gl/iI3Y0. Table 2 demonstrates the running time efficiency of computing 
/gp,M for random polynomials / and random integers M chosen as before but for 
relatively large n and 2d and with sparsity conditions on the size of 

Table 1. Average running time for f gPt M and f^l M (seconds) 



n 


2d 


4 6 8 10 


3 


/gp,M 

AO) 

Jsos,M 


0.03 0.09 0.96 4.73 
0.05 0.56 6.42 62.28 


4 


/gp,M 

AO) 

J sos. A/ 


0.04 0.89 34.90 278.43 
0.16 7.74 154.17 


5 


/gp,M 

AO) 

J sos,M 


0.10 8.25 48.28 1825.56 
0.53 69.49 



Table 2. Average running time for f gp ,M (seconds) for various 
constraints on 



n 


2d\\n(f)\ 


10 


20 


30 


40 


50 




20 


0.52 


0.62 


1.91 


4.36 


5.63 


10 


40 


0.75 


1.42 


2.1 


5.08 


11.16 




60 


0.86 


1.72 


3.1 


6.48 


13.07 


20 


20 


3.69 


18.11 


17.11 


44.78 


46.51 


40 


3.75 


18.82 


37.52 


59.55 


114.05 




60 


7.31 


27.33 


46.05 


96.86 


164.56 


30 


20 


3.16 


19.63 


34.81 


44.04 


175.5 


40 


6.07 


22.72 


105.77 


217.07 


315.85 




60 


13.71 


72.81 


132.04 


453.05 


667.87 


40 


20 


6.67 


37.22 


63.09 


131.03 


481.71 


40 


11.21 


76.03 


83.91 


458.75 


504.6 




60 


24.97 


114.45 


355.56 


796.52 


1340.76 



hardware and Software specifications. Processor: Intel® Core™2 Duo CPU P8400 @ 
2.26GHz, Memory: 2 GB, OS: Ubuntu 12.04-32 bit, Sage-4.8 
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We compare values of / gp ,M with corresponding values of /^J M for various 
choices of / and M. 

Example 4.1. Let / = w 6 + x 6 + y 6 + z 6 + 7w 4 y - I0w 3 xy + 5wx 3 y - 3w 3 y 2 - 
3w 2 xy 2 + 9wxy 3 - 10xy 4 + 7w 4 z + wx 3 z - bxyz 3 - 5z 5 + 8w 4 + 8w 2 x 2 - 4wx 3 - 
w 3 y + 2wx 2 y + 3w 2 y 2 - wxy 2 + wy 3 + 7w 2 xz - 3y 3 z + w 2 z 2 + 2y 2 z 2 - 2w 3 + 8x 3 - 
5w 2 y + 8x 2 z + 3xz - 3z + 5, then: 

/ gP ,i » -39.022 £l A » -5.519 

/gp.io « -213.631 ifi^ » -67.947 

/ gp , 10 2 « -1215.730 /^^ « -489.009 

/ gp w -9580211.794 / sos 'w -458107.262 

Example 4.2. Let / = 8w 6 +6x e +4y e +2z 6 -3w 3 x 2 +8w 2 xyz-9xz 4 +2w 2 xz-3xz 2 1 
then 

/ gP ,i « -6.605 
/g P ,io « -27.151 
/ gP ,io 2 « -73.458 
« -74.971 

Example 4.3. For / = -7x 3 y 4 + 13x 2 y 5 + 5y 4 z + 18xz 4 - 5z 2 with 2d = 8 

« -23.4559 / s ( o s } ;1 « -19.4797 
/gp.io « -117.9727 10 « -92.6547 

/gp.io^ « -736.0259 /J^ « -668.221 

We can compute / gPl M in cases where computation of M breaks down. 

Example 4.4. For / = -9w 12 a; 9 y 12 z 5 + 19w 8 a; 2 yz 20 - 3w n a; 6 y 9 z 4 - 3w 13 x 14 z - 
18w 4 x 12 y 3 with 2d = 40 

/ gp ,i « -20.0645 
/gp.io « -106.4946 
/ gPil0 2 « -584.027 

Example 4.5. For 

/ = YhLo X T + X%X%X^X 7 x\x 9 X W x\ 1 Xl2 - Ux^XzX^X-jxlxiQX^X^XifiXisXig + 
19xlxfxlx 9 Xi2X 2 7 Xi8xl 9 -10XoXiX2xlxi2Xi 5 Xi 7 xl s X 4 9 -llx 2 t X2X 3 X 5 XQX 4 2 X 4 5 Xi6Xi7 + 

15xIxIxqX$x 9 x1 4 x 4 5 x1 8 x 2 9 + 2xix 2 l x 3 x e xi xl 1 xi 3 xi 5 xi 7 xi$xl g , 

/gp.io « -41.6538 
/ g p,io= « -340.6339 
/ gP ao3 « -2774.217 
/ gp w -84853211002.07141 



/Sio » -27-151 
/ sos « -74.971 
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5. Proof of Theorem 2.4 

With A > fixed, let us apply Proposition 2.1 to the polynomial fx € M[x] 2t i, so 
as to obtain a lower bound G(A) on G(A) defined in (5). Then G(A) :— f\{0) — p\, 
with 



p\ = min (2d — 

a£A(f x )<™ 



2d \z, 



l/(2d-|«|) 



s.t. 22 z a,t < (hhd,i, i = l,...,n 

«£A(/ X ) 



2,1 



(mt) (^r = 1 ! a€A(A),N=24 
Notice that f2(/^) = and (f\) a = f a for all a £ fi(/). Moreover, 

/a(0) = /(0) - AM; (/ A ) 2<M = f 2d>i + A, Vi = 1, . . . ,n. 
And so, with A > 0, G(A) := /(0) - MA - p A , with 



Pa = min ^ (2d-|a|) 

q£A(/)< m 



M M / x al V(2rf-I«l) 

Jo \ / a 
2d 



s.t. z a>i < / 2 d,i + A, i = 1, . . . , n 

«6A(/) 

(f) (l?) =1 ! «eA(/), |a|=2d, 
is a lower bound on G(A) for every A > 0. Next, recall that 

J2d,l > J2d,2 > ••• > J2d,n- 

Let 

A := max{0, -f 2 d,n}- 

For < A < A , f 2d , n + A < 0, so p A = oo, i.e., G(A) = -co. For A > A , Pa < P\ , 
i.e., G(A) > G(A ) - M (A - A ). Consequently, 

maxG(A) = max G(A) = maxG(A). 

A>0 A>A A>A 

For A > Ao, using the new variables := f 2dl i + A > 0, i = 1, . . . , n, one has: 
Ui = Ui-i - {f 2 d,i-i - hd,i), i = 2,...,n, 

or equivalently, 



Ui f 2d ,i-i — f: 



2d., 



Ui-1 



Ui-l 



1, i = 2, . . . , n. 
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In addition, the constraints z a ^ < fed.i + A, read 



Ui 



V — < 1, i = l,...,n. 

£ ^ 7/ . 

Finally, as A = u\ - f 2 d,i, then /(0) - MX = /(0) + M/ 2d ,i - Mu x . Therefore, for 
A > Ao fixed, and 

\ . r u i , f 2 d,i-X — J 2d,: 

ui = \ + f 2( t,i; 1 

G(A) = /(0) + M/ad,! - M (u) with 



= 1, i = 2,. . . ,n, 



~>M 



(u) = Mu x + min ^ {2d - \a\) 

aeA(f)< 2d 



2d 



V(2d-|a|) 



s.t. < 1, i = l,...,n 

aeA(/) " 4 

(l!) 2 " CST = !; «GA(/), |a|=2d. 



And so, 



where 



maxG(A) = max G(A) = maxG(A) = /(0) +Mf 2d ,i - Pm, 

A>0 A>Aq A>Aq 



= nun Mui + > (2d - 

z a .u ^ — ' 

aEA(/)< 2i 



2(/ 



l/(2d-|a|) 



S.t. 



(9) 



I /2d,i-l~/2d,z 
lii 1 Ui — 1 

Ml — 



1, i = 2, . . . , n 



E 



< 1, i — 1, . . . , n 



aeA(f) 
/ \ 2d 



if)" = 1; a E A(f), |o| = 2A 

Notice that (9) is not a geometric program because of the presence of n — 1 
posynomial equality constraints. To obtain a geometric program, observe that in 
(9) we can relax the n — 1 posynomial equality constraints 

Ui f2d,i-l — f2d,i 



(10) 

to the posynomial inequality constraints 

Ui f2d,i-X — f2d,; 



(11) 



Ui-x 



Ui-x 



= 1, i = 2,. 



< 1, i = 2,...,n, 
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without changing the optimal value. Indeed, suppose that u is an optimal solution 
of (9) with (11) in lieu of (10). Then increase u-i to u' 2 := U2 + 8 2 with 5 2 > so 
that 

U 2 + S 2 f<2d,\ — /2d,2 _ ^ 



Ui Ui 

Since < u 2 < uL the constraint ^ < 1 and the constraint + hd - 2 ~ hd - 3 < 

£ — ^ yf u 2 u 2 

aeA 2 

1, are satisfied. Therefore, one may repeat the process now with M3, i.e., increase 
W3 to w 3 = U3 + 83 with <53 so that 

U3 + ^3 /2d,2 — hdfi _ , 

u' 2 u' 2 

Since < u 3 < uL the constraint < 1 and the constraint /2d ' 3 ~; /2d - 4 < 

1, are satisfied, etc. Iterate the process to finally obtain a feasible solution ((z a ,i), u') 
for (9), with the desired property. In addition, since wi and (z a ,i) have not been 
changed, the cost associated to the new feasible solution ((z a ,i), u') is the same. □ 
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